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In the framework of a Fermi liquid theory it is considered the possibility of ferro- 
magnetic and antiferromagnetic phase transitions in symmetric nuclear matter with 

cd ' 

Skyrme effective interaction. The zero temperature dependence of ferromagnetic 
and antiferromagnetic spin polarization parameters as functions of density is found 
for SkM*, SGII effective forces. It is shown that in the density domain, where both 
type of solutions of self-consistent equations exist, ferromagnetic spin state is more 
preferable than antiferromagnetic one. 
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Q 1 I. INTRODUCTION 

a 

The appearance of spontaneous spin polarized states in nuclear matter is the topic of 
a great current interest due to relevance in astrophysics. According as nuclear matter is 
spin polarized or not, drastically different scenarios of supernova explosion and cooling of 
neutron stars can be realized. The possibility of a phase transition of normal nuclear matter 
to ferromagnetic state was studied by many authors. In Ref. [1] within a hard sphere 
gas model it was shown that neutron matter becomes ferromagnetic at p m 0.41 fm -3 . In 
Refs. was found that the inclusion of long-range attraction significantly increases 

the ferromagnetic transition density (e.g., up to p ~ 2.3 fm~ 3 in the Brueckner theory with 
a simple central potential and hard core only for singlet spin states [3]). Calculations of 
magnetic susceptibility in Ref. |4j with Skyrme effective forces showed that a ferromagnetic 
transition occurs at p ~ 0.18-0.26 fm~ 3 . The Fermi liquid criterion for the ferromagnetic 
instability in neutron matter with Skyrme interaction is reached at p ~ 2-4^ [5], Po being 
nuclear matter saturation density. In Ref. [fj there were formulated general conditions on the 
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parameters of neutron-neutron interaction, which result in the magnetically ordered state 
of neutron matter. The spin correlations in dense neutron matter within the relativistic 
Dirac-Hartree-Fock approach with the effective nucleon-meson Lagrangian were studied in 
Ref. pi, predicting the ferromagnetic transition at several times nuclear matter saturation 
density. The importance of the Fock exchange term in the relativistic mean-field approach 
for occurrence of ferromagnetism in nuclear matter was established in Ref. . The stability 
of strongly asymmetric nuclear matter with respect to spin fluctuations was investigated 
in Ref. 9j, where it was shown that even small proton admixture favors the ferromagnetic 
instability of the system. This conclusion was confirmed also by calculations within the 
relativistic Dirac-Hartree-Fock approach to strongly asymmetric nuclear matter [lo| . 

If to consider the models with realistic nucleon-nucleon (NN) interaction, the ferromag- 
netic phase transition seems to be suppressed up to densities well above g lj] I3- In 
particular, no evidence of ferromagnetic instability has been found in recent studies of neu- 
tron matter and asymmetric nuclear matter within the Brueckner-Hartree-Fock 
approximation with realistic Nijmegen II, Reid93 and Nijmegen NSC97e NN interactions. 
The same conclusion was obtained in Ref. Q], where magnetic susceptibility of neutron 
matter was calculated with the use of the Argonne v is two-body potential plus Urbana IX 
three-body potential. 

Here we will continue the study of spin polarizability of nuclear matter with the use 
of effective NN interaction. As a framework of consideration, it is chosen a Fermi liquid 



(FL) description of nuclear matter 13> UM- As a potential of NN interaction, we use the 



Skyrme effective interaction, utilized earlier in a number of contexts for nuclear matter cal- 
culations Since calculations of magnetic susceptibility with the Skyrme effective 
forces show, that nuclear matter undergoes the ferromagnetic phase transition at some criti- 
cal density, it is quite natural step to determine the density dependence of the ferromagnetic 
spin polarization parameter of nuclear matter. Besides, in this study we explore the pos- 
sibility of antiferomagnetic phase transition in nuclear matter, when the spins of protons 
and neutrons are aligned in the opposite direction. The antiferromagnetic spin polarization 
parameter will be determined as a function of density as well. Then we study the question 
of thermodynamic stability of ferromagnetic and antiferromagnetic spin states and clarify 
which phase is thermodynamically preferable in the density region, where both solutions of 
self-consistent equations exist simultaneously. 
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Note that we consider thermodynamic properties of spin polarized states in nuclear matter 
up to the high density region relevant for astrophysics. Nevertheless, we use a pure nucleonic 
description of nuclear matter, although other degrees of freedom, such as pions, hyperons, 
kaons, or quarks could be present at such high densities. 



II. BASIC EQUATIONS 

Normal states of nuclear matter are described by the normal distribution function of 
nucleons f KlK2 = Tr ga+ 2 a Kl (k = (p,cr, r), p is momentum, er(r) is the projection of spin 
(isospin) on the third axis, g is the density matrix of the system). The energy of the system 
is specified as a functional of the distribution function f,E = E(f), and determines the 
single particle energy 

««-(/) = f^- (i) 

The self-consistent matrix equation for determining the distribution function / follows from 



the minimum condition of the thermodynamic potential 



17| and is 



/ = {exp(y £ + n) + 1}^ = {exp(F O + I} -1 • (2) 

Here the quantities e,Y 4 are matrices in the space of k variables, with Y^ KlK2 = Y^ T1 5 niK2 
(t\ = p,n), Y = 1/T, Y ip = —fip/T and Y in = —fi^/T the Lagrange multipliers, /i° 
and /x° the chemical potentials of protons and neutrons and T the temperature. We shall 
study the possibility of formation of different types of spin ordering (ferromagnetic and 
antiferromagnetic) in nuclear matter. 

The normal distribution function can be expanded in the Pauli matrices a, and in spin 
and isospin spaces 

f(p) = /oo(p)o"oT + /3o(p)^3r (3) 

+ /03(p)o-07- 3 + /33(P)CT3^3- 

For the energy functional, being invariant with respect to rotations in spin and isospin spaces, 
the structure of the single particle energy is similar to that of the distribution function /: 

e(p) = £oo(p)o-o7o + e 30 (p)a 3 T (4) 
+ € 03 (p)croT 3 + e 33 (p)a 3 T 3 . 
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Using Eqs. (J2J), Q, it is possible to express evidently the distribution functions 
Zoo, /so, /03, /as in terms of the quantities e: 

/oo = ~{n(u; +)+ ) + 7i(w + _) + 7i(w_,+) + n(w_,_)}, 
/so = ^{n(u; +)+ ) + _) - n(w_ +) - 
/os = J{n(w + ,+) - n(u +t J) + - n(w_,_)}, 

/as = ~{n(a; + , + ) - n(w +> _) - n(w_ +) + n(w_,_)}. (5) 
Here n(a;) = {exp(F ^) + I} -1 and 

^+,+ — £00 + £30 + £03 + £33, 

UJ +- — ^OO + ^30 — £o3 — £33, 
^-,+ = ^OO — £30 + £o3 — ^33, 
U-,- = ^00 — ^30 — 6)3 + ^33, 



where 



^00 — ^00 — A^OO) £30 _ £ 30, 

C - c 

"503 — £ 03 ~ ^035 S33 — ^33; 

,0 1 ,,0 ,,0 ,,0 



,0 



/% ~~ o ' ^03 



2 2 

As follows from the structure of the distribution functions /, the quantity oo± t ±, being the 
exponent in Fermi distribution function n, plays the role of the quasiparticle spectrum. In 
the considering case the spectrum is four-fold split due to spin and isospin dependence of 
the single particle energy e(p) in Eq. (jlj). The distribution functions / should satisfy the 
normalization conditions 

4 x - 

y2^/oo(p) = £, ( 6 ) 
p 

4 x - 

y 2^ ^ 03 ( p ) = Qv~ Qn = -ag, (7) 
p 

4 x - 

y 2^ /so(p) = &[-Qi = A^ TT , (8) 
p 

y /3s(p) = (ft* + eni) - + 0nt) = A £u- ( 9 ) 

p 
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Here a is the isospin asymmetry parameter, g p ^, g p ^ and g n ^, g n ^ are the proton and neutron 
number densities with spin up and spin down, respectively; g^ = g P i + g n ^ and g± = g P i + g n ± 
are the nucleon densities with spin up and spin down. The quantities Ag^ and Ag^ may be 
regarded as ferromagnetic (FM) and antiferromagnetic (AFM) spin order parameters: if all 
nucleon spins are aligned in one direction (totally polarized FM spin state), then Ag^ = g 
and Apjj = 0; if spins of all protons are aligned in one direction and spins of all neutrons 
in the opposite one (totally polarized AFM spin state), then Ag^ = g and Ag^ = 0. 

To obtain the self-consistent equations, it is necessary to specify the energy functional of 
the system, which we write in the form 

E(f) = E (f) + E mt (f), (10) 

£ (/) = 45> (p)/oo(p), £ o(p) = £-> 

p 

Eint(f) = 2^{£ 0(P)/00(P) +^30(P)/30(P) 
P 

+ ^03(P)/0S(P) + £33(P)/33(P)}, 



^oo(p) = C/o(k)/oo(q), k = 

q 

^o(p) = ^$>i(k)/ 3 o(q), 
q 

£o3(p) = ^^C/ 2 (k)/o3(q), 
q 

^33(p) = ^E f/ 3(k)/ 33 (q). 
q 

Here mo is the bare mass of a nucleon, Uq(]s), L^k) are the normal FL amplitudes, 
^oo, ^30, ^03, ^33 are the FL corrections to the free single particle spectrum. With allowance 
for Eqs. ((TJ) and (fTUjl. we obtain the self-consistent equations in the form 

£oo(p) = e (p) + e 00 (p) - Ho , £ 30 (p) = £3o(p), (11) 

£o 3 (p) = £03(p)-/^03> ^3 3 (P) =£ 3 3(P)- 

Further for obtaining numerical results we shall use the Skyrme effective interaction. In the 



case of Skyrme forces the normal FL amplitudes read 



U (k) = 6t + + ^[3h + £ 2 (5 + 4x 2 )]k 2 , (12) 
C/i(k) = -2t (l - 2x ) - ^3^(1 - 2x 3 ) 

- -2^(1 - 2x x ) - t 2 (l + 2x 2 )]k 2 = a + 6k 2 , 



ft 

C/ 2 (k) = -2t (l + 2a;o) - -t 3 ^(l + 2a; 3 



f [t 1 (l + 2a; 1 )-t 2 (l + 2x 2 )]k 2 , 



1. « 2 



C/ 3 (k) = -2t - i* 3 / - 72 (*i - t 2 )k 2 = c + rfk 2 , 

o ft 

where ti,Xi,/3 are phenomenological constants, characterizing the given parametrization of 



23( and SGII |2J potentials, 



Skyrme forces. In numerical calculations we shall use the SkM* 
developed to fit the properties of systems with small isospin asymmetry. With account of 
the evident form of FL amplitudes and Eqs. (JHJ)-©, one can obtain 

£oo = tt~ A*00) (13) 

2m 00 

p 2 

e 3 o = (a + 4)^f + 4 (q2)30 ' (15) 

fe = ( C +4)^r + l (q2)33 ' (16) 

where the effective nucleon mass moo and effective isovector mass mo3 are defined by the 
formulae: 

h 2 h 2 g 



2m o 2m 16 

,2 



+ ^[3ti+t 2 (5 + 4x 2 )], (17) 



n ' f? [t 1 (l + 2x 1 )-t 2 (l + 2x 2 )], 



2m 03 16 

and the renormalized chemical potentials /%)>A*03 should be determined from Eqs. ©, 
(J7J). In Eqs. (|T3|). (fTB^) (q 2 )3o, (q 2 )33 are the second order moments of the corresponding 
distribution functions 

(q 2 >3o = ^£q 2 / 3 o(q), (is) 

q 

<q 2 > 33 = ^;>> 2 /33(q)- (19) 

q 
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Thus, with account of the expressions (0) for the distribution functions /, we obtain the 
self-consistent equations (fTSj). (jTUj) for the effective chemical potentials /ioo,/-*03, 

FM 

and AFM spin order parameters A^^, Ag^, and the second order moments (q 2 )3o, (o. 2 )33- 

III. FERROMAGNETIC AND ANTIFERROMAGNETIC SPIN ORDER 
PARAMETERS AT ZERO TEMPERATURE 

The early researches on spin polarizability with Skyrme effective interaction were based on 
the calculation of magnetic susceptibility and finding its pole structure 0,0], determining the 
onset of instability with respect to spin fluctuations. Here we provide the direct calculation 
of FM spin polarization as a function of nuclear matter density at zero temperature. Besides, 
we study the possibility of AFM spin ordering in nuclear matter and competition between 
these two types of ordering. 

Let us consider the zero temperature behavior of spin polarization in symmetric nuclear 
matter (g p = g n ). The FM spin ordering corresponds to the case Ag^ ^ 0, (q 2 )3o ^ 
0, Agfj = 0, (q 2 )33 = 0, while the AFM spin ordering to the case ^ 0, (q 2 )33 ^ 

0, Agf-f = 0, (q 2 )3o = 0. In the totally ferromagnetically polarized state nontrivial solutions 
of the self-consistent equations have the form 



Here kp = (37i 2 g) 1 ^ 3 is Fermi momentum of symmetric nuclear matter in the case when 
degrees of freedom, corresponding to spin up of nucleons, are open while those related to 
spin down are inaccessible. For totally antiferromagnetically polarized nuclear matter we 
have 



The Fermi momentum A;^ is given by the same expression as in Eq. (ffllj) since now de- 
grees of freedom, related to spin down of protons and spin up of neutrons, are inaccessible. 
The results of numerical determination of FM Ag^/g and AFM Ag^Jg spin polarization 
parameters are shown in Fig. ^ for the SkM* and SGII effective forces. 

The FM spin order parameter arises at density g ~ 2g Q for the SkM* potential and at 
q ~ 2.75^o for the SGII potential. The AFM order parameter originates at g ~ 3.3^o for the 
SkM* force and at g ~ 5go for the SGII force. In both cases FM ordering appears earlier 



Ag n = g, (q 2 ) 3 o 




(20) 



Aqu = q, (q 2 )33 




(21) 
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than AFM one. Nuclear matter becomes totally ferromagnetically polarized (Ag^/g = 1) 
at density g ~ 2.7^o for the SkM* force and at g ~ 3.9^o for the SGII force. Totally 
antiferromagnetically polarized state (Ag^Jg = 1) is formed at g ~ 4.5^o fo r the SkM* 
potential and at g ~ 7.2^ for the SGII potential. 

Note that the second order moments (q 2 )3o, (q 2 )33 of the distribution functions /30, ^33 
play the role of the order parameters as well. In Fig. |2] it is shown behavior of these 
quantities normalized to their value in totally polarized state. The ratios 5(q 2 )3o/3^»fc 2 ? and 
5 (q 2 ) 33/3gkp are regarded as FM and AFM order parameters, respectively. The behavior of 
these quantities is similar to that of the spin polarization parameters in Fig. with the same 
values of the threshold densities for appearance and saturation of the order parameters. 

In the density domain, where FM and AFM solutions of self-consistent equations exist 
simultaneously, it is necessary to clarify, which solution is thermodynamically preferable. 
With this purpose it is necessary to compare the free energies of both states. The results 
of the numerical calculation of the free energy density, measured from that of the normal 
state, are shown in Fig. El One can see, that for all relevant densities FM spin ordering is 
more preferable than AFM one, and, moreover, the difference between corresponding free 
energies becomes larger with increasing density, so that there is no evidence, that AFM spin 
ordering might become preferable at larger densities. 

In conclusion, we have considered the possibility of spontaneous appearance of spin po- 
larized states in symmetric nuclear matter, corresponding to FM and AFM spin ordering. 
The study has been done in the framework of a Fermi liquid description of nuclear mat- 
ter, when nucleons interact via Skyrme effective forces (SkM*, SGII potentials). Unlike the 
previous considerations, where the possibility of formation of FM spin polarized states was 
studied on the base of calculation of magnetic susceptibility, we obtain the self-consistent 
equations for the FM and AFM spin polarization parameters and solve them in the case of 
zero temperature. It is shown that FM order parameter appears at densities 2-2.75^o and 
AFM one at densities 3.3-5^o- in the density region, where both type of solutions exist, FM 
spin ordering wins competition for thermodynamic stability. 

The author acknowledges the financial support of STCU (grant No. 1480). 
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FIG. 1: FM and AFM spin polarization parameters as functions of density at zero temperature for 
(a) SkM* force and (b) SGII force. 
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FIG. 2: The same as in Fig. ^ but for the second order moments of the distribution functions. 
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FIG. 3: Free energy density, measured from that of the normal state, for FM and AFM spin 
ordering as function of density at zero temperature for (a) SkM* force and (b) SGII force. 



